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Abstract 

In July 1967 the Explorer 35 spacecraft was launched into a lunar 
orbit to study the interaction of the magnetized solar wind plasma with 
the moon. The experimental results indicate that a detached bow shock 
wave does not exist in the vicinity of the moon. Thus the flow condi- 
tions near the moon do not resemble those near the earth's magnetosphere 
In this paper the solar wind flow around the moon is treated theoreti- 
cally as a free molecule flow of magnetized plasma. Analytical results 
are obtained to describe the distribution of the ion density and the 
ion flux in the vicinity of the moon. A sizable empty cavity is pro- 
duced in the near-downstream vicinity of the moon; across the boundary 
of the cavity the plasma flow changes rapidly from the undisturbed con- 
dition to the void condition. The disturbed region in the downstream 
forms a long wake on the dark side of the moon. The direction of the 
wake is deflected slightly toward the field lines. The thickness of 
the wake, measured perpendicular to the plane of the solar wind velocity 
and the interplanetary magnetic field, is constant at one lunar diameter 
On the other hand, the width of the wake, measured parallel to that 
plane, increases with distance from the moon. The length of the wake 
is of the order of 100 lunar radii. 



1. Introduction 


The solar wind ^ 2 flows approximately radially outward 
from the sun in all directions continuously. It flows past the 
earth’s orbit at a supersonic velocity. The solar magnetic lines 
of force are carried outward by the solar wind and are twisted into 
the form of Archimedean spirals. The earth’s magnetosphere presents 
a blunt-nosed obstacle in the course of the solar wind blocking the 
motion of both the solar wind plasma and the interplanetary magnetic 
field lines. In the vicinity of the earth's magnetosphere, the 
solar wind plasma behaves as a continuum, and a detached bow shock 
wave is formed standing on the sunward side of the earth's magneto- 
sphere. Satellite measurements 3 ' 4 have detected the presence of 
the detached bow shock and the magnetosheath region between the mag- 
netosphere and the bow shock. Theoretical calculations 5-7 have been 
successful in explaining the gross features and the overall shape of 
the interaction region. As the supersonic solar wind interacts with 
the moon, it might be expected that the flow conditions near the 
moon would resemble those near the earth's magnetosphere with a de- 
tached bow shock standing on the sunward side of the moon. On July 
22, 1967 the NASA Explorer 35 spacecraft was placed into a lunar 
orbit to study the interaction of the solar wind with the moon, but 
no detached bow shock has been observed 8 ' 9 . The experimental results 
indicate that the actual interaction of the solar wind with the moon 
is quite different from that of the solar wind with the earth. 

In ordinary fluid mechanics, a division of fluid flows into 
various regimes is determined by the ratio of the collisional mean- 
free-path X to the characteristic length of the flow L. In the 
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continuum flow regime, where \ « L, particle-particle interactions 
dominate over particle- surface interactions. In the free molecule 
flow regime, where X » L, particle-surface interactions dominate 
over particle-particle interactions. Between these two limiting 
regimes is the transition flow regime. When a supersonic flow en- 
counters a blunt-nosed obstacle, in the continuum flow regime a 
detached shock wave is generated upstream of the obstacle, while in 
the free molecule flow regime no shock wave can be formed in the 
vicinity of the obstacle. Absence of a detached bow shock in the 
vicinity of the moon indicates that the interaction of the solar 
wind plasma with the moon behaves like a free molecule flow rather 
than a continuum flow. 

The collisional mean-free-path of the solar wind plasma 
near the earth's orbit is of the order of 1 AU. In lieu of inter- 
particle collisions, the interplanetary magnetic field has the 
important effect of causing the solar wind plasma to behave as a 
continuum for the interaction between the solar wind and the earth's 
magnetosphere. Thus the interaction length must be small in compar- 
ison with the dimension of the earth's magnetosphere. The exact 
scale of this interaction length is not well understood yet. The 
diameter of the earth's magnetosphere is about 50 times that of the 
moon, which suggests that the small size of the obstacle is responsible 
for the absence of continuum phenomena near the moon. In addition 
to the difference in their sizes, their contrary responses to the 
impinging plasma and field lines also affect the flow patterns around 
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the moon and the earth's magnetosphere. Firstly, the blunt nose of 
the earth's magnetosphere is very effective in blocking the pene- 
tration of the interplanetary magnetic field lines. On the contrary, 
the experimental results 8 obtained from Explorer 35 indicate that 
the moon does not have an intrinsic magnetic field and the average 
electrical conductivity of the moon is very low, the impinging field 
lines appear to diffuse through the lunar body very rapidly. Secondly, 
all the charged particles striking on the moon are captured by the 
cold surface of the moon; they will be neutralized and then re-emitted 
as cold neutral particles. Whereas the blunt nose of the earth's 
magnetosphere bounces off all impinging solar wind particles. 

In this paper, we will treat the motion of ion particles 
of the solar wind in the vicinity of the moon as a free molecule 
flow of magnetized plasma. The moon's surface will be considered as 
perfectly absorbing with respect to the impinging charged particles. 

2 . Ana lysis 

Let us consider that in the undisturbed upstream flow the 
solar wind moves with a steady uniform velocity U Q with respect to 
the center of the moon, and a steady uniform interplanetary magnetic 
field B 0 is carried with the solar wind (Fig. 1). The velocity of 
a charged particle in the upstream can be expressed as the vector 
sum of a guiding-center velocity Vq and a Larmor motion velocity C ± . 

The component of the guiding-center velocity normal to the field 
lines is the electric drift, which is exactly the component of JJp 
normal to B 0 . All charged particles have the same electric drift 
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velocity, which means that the guiding-centers do not have thermal 
motion transverse to the field lines. The component of the guiding- 
center velocity parallel to the field can be expressed as the sum 
of a thermal velocity J3|j and a mean mass flow velocity which is 
equal to the component of parallel to B 0 . Thus we can express 
the guiding center velocity as composed of two parts: the mean 

mass flow velocity y o and a thermal motion velocity along the field 
lines J2„ (Fig. 2): 

% = JIo + C,, (1) 

Physically, this equation means that the guiding center of a charged 
particle is constrained to move along a magnetic line which is itself 
moving with a velocity U Q . 

The ion particles are rapidly gyrating about the guiding 
center with a Lamor radius of the order of 80 km, which is only a 
small fraction of the moon's radius (Elm = 1738 km). In order to 
find the lowest order solution for the ion flow in the vicinity of 
the moon we approximate each ion particle by an artificial guiding* 
center particle a guiding-center particle is considered as 

a charged particle located at its guiding center, moving with its 
guiding-center velocity, and with the oscillation about the guiding 
center representing internal energy, magnetic moment, and angular 
momentum of the particle. A collection of artificial guiding-center 
particles behaves as a one -dimensional gas which has thermal motion 
only along the field lines. 
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Next we assume that the distribution function of this one - 
dimensional guiding-center gas is Maxwellian in the undisturbed 
upstream flow, namely, 

f G = n O “|| ex P (-Of | cj ) (2) 

Here ajj = m^/ (2k T )t ), with the parallel temperature Tj| 
defined by 

k T„ 12 = (l/n 0 ) J.“ (mcj /2)f G dC„ 

In the free molecule flow in the vicinity of the moon, the 
particle-surface interaction dominates over the particle-particle 
interaction. We make the additional assumption that when the solar 
wind plasma arrives in the vicinity of the moon, each guiding center 
particle continues to move with its initial guiding-center velocity 
along a straight-line path unless it is captured by the moon’s cold 
surface. The guiding-center distribution function remains constant 
on each straight path. Inevitably, the magnetic and the electric 
fields will be slightly disturbed in the moon’s wake region. The 
drifts due to non-uniformity of the disturbed fields are small com- 
pared with the initial guiding-center velocity, and these "first- 
order" drifts are negligible in a "zero-order" solution. 

Based on the assumptions outlined above, we can now calculate 
the ion density and the ion flux in the vicinity of the moon. We 
will choose a coordinate system moving with the center of the moon, 
with its origin located at the moon’s center, the X-axis parallel 
to the velocity vector U 0 , and the Z-axis along the direction of 
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Bo x Uo (Fig. 1). The scale of the coordinate system is normalized 
by using the moon's radius as a unit length. 0 denotes the direction 
angle of the interplanetary magnetic field. e x and e^ denote the 
unit vectors along the X-direction and the B 0 -direction respectively. 

The XY-plane which is parallel to both the solar wind velocity 
U 0 and the interplanetary magnetic field B c is a plane of symmetry 
for the flow around the moon. The guiding-center velocity has no 
component along the Z-direction. Thus the motion of each guiding - 
center particle is restricted on a Z^onstant plane. The moon’s 
surface would not intercept any guiding-center particle's trajectory 
in the region |z| >1, the flow field is undisturbed in this region. 

In the region |z| £ 1, we may first draw a line on the Z= constant 
plane tangential to the moon's surface on its upstream side (Fig. 3). 
At any point on the upstream side of this line, we can find guiding- 
center particles distributed over the whole spectrum of thermal 
velocity (i.e., -« < < + <»)> the flow field is not disturbed due 

to the presence of the moon. The ion density is n Q and the ion flux 
is n 0 U 0 in the undisturbed region. The flow in the downstream side 
of the first tangential line is more or less disturbed. 

If we draw a second line on the Z=constant plane tangential 
to the moon on its downstream side, we can divide the disturbed flow 
into the three regions identified as region I, II and III in Fig. 3. 

At any point P(r) in the disturbed region, the ion density and the 
ion flux can be calculated if the guiding center distribution 
function is known at that point as a function of C|| . When an un- 
intercepted particle arrives at P, it still moves with its initial 
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JJ 0 and , and carries the same Maxwelliam distribution function 
f g with it* Because the moon’s surface behaves as a perfectly 
absorbing surface with respect to charged particles, no particle 
can arrive at point P if its straight-line trajectory has been in- 
tercepted by the moon’s surface. Therefore at point P the guiding- 
center distribution function is Maxwellian in a certain interval 
in the one-dimensional -space, and is zero outside of this 
interval* Let Ij denote the interval of in region j within 
which the distribution function is Maxwellian, We can express 
these intervals as 

Ii: (- ro > - Y1 u o> 

x 2 : (V2 u o> + 00 ) ► (3) 

1 3 : (- ^ ** Y\ U 0 ) and (y 2 u o> 00 ) 
where Y1 88 sin (o' + (3)/ sin (0 - o' - 0) 

= sin (a - (3)/ sin (0 + or - 0) 
ot = arcsin ((1 - Z^)^(X^ + Y^) 

0 = arctan (Y/X) 

Now the guiding-center distribution function in the disturbed 
regions (regions I, II and III) can be expressed as 

f G = TT’* n c Q? (| exp (-c*2 Cj* ) for C,, e Ij 
= 0 for C i Ij 

The ion density, n 0 , and the ion flux, nu, can be calculated from 



n " J I j f G dC ll 

n H. " J* (U 0 + £„ ) f G dc || 
j 


and 


( 6 ) 
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Let S = U/(2kT n /m^ denote the so-called "speed ratio" in the 
free molecule flow regime. The integration, from (5) and (6), 
yields the following forms: 

In region I, 

n/n Q = (%) erfc(Y]S) 

nu/(n 0 U 0 ) = (%) [e x erfc CY]S) - e^ exp (-y 1 2 S 2 ) / (n^S) ] 

In region II, 

n/n Q =(1/2) erfc (Y 2 S) 

nu/(n 0 U 0 ) = (%) [e x erfc (Y 2 S) + e x exp (-Y^S 2 ) /(tt^S) ] 

In region III, 


1 ( 8 ) 


n/n Q = (1.) [erfc (YjS) + erfc (Y 2 S )1 
nu/(n 0 U 0 ) = (%) e x [erfc (Y]S) + erfc (Y 2 S)] - 
ei [exp ( - Yl 2 S 2 ) " ex P (~Y 2 2 & ]/@n^S) 


V (9) 


Erfc(^) is the complementary error function of ^ defined as 


erfc(^) = 1 - erf (|) 




dt 


3. Results 

At the earth's orbit, the typical interplanetary conditions 
of the past few years directly observed from space satellites are 
of the order of U - 400 km/sec, T — 10“* °K, n ~ 5 ions/cm, B ~ 6y 
and 0 ~ 135° or 315°. The temperature parallel to the magnetic 

field, Tj| , is about twice the temperature transverse to the mag- 
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netic field, T x . Based on these values some numerical computations 
have been carried out. 
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We shall first discuss the ion flux distribution. Fig. 4 shows 
the distribution of ion flux on the XY-coordinate plane for S=5 and 
0 = 90°, 105°, 120°, 135° and 150°. The guiding-center particles are 
insensitive to the polarity of the field line; when 0=0^ + 180° the 
flow pattern is exactly the same as that with 0 = 0 -^. The ion flux 
is disturbed mainly downstream from the moon. A sizable empty cavity 
is produced in the near -downstream vicinity of the moon; across the 
constant flux contours around the cavity the ion flux changes rapidly 
from the undisturbed condition to the void condition. A long wake is 
formed on the dark side of the moon. These features agree with the 
observations reported in references 9 and 14. When the field lines are 
perpendicular to the direction of the relative velocity of the solar 
wind, the wake region is relatively short; the length of the wake re- 
gion increases as the acute angle between the two directions decreases. 
As 0 tends to 180°, the shape of the empty cavity approaches a circular 
cylinder with a radius of 1 R m ; the ion flux and the ion density are 
undisturbed outside the cylinder and zero inside. Fig. 5 plots the 
distribution of ion flux on the XY-coordinates plane for 0 = 135° and 
S = 3, 7 and 10. It is seen that the moon's wake becomes narrower and 
longer as the speed ratio S increases. 

Let 9 denote the angle between the direction of the ion flux and 
that of the solar wind velocity with respect to the moon, 

0 = arctan (uy/u x ) . 

The constant 9 contours are plotted in Fig. 6 to show the change of ion 
flux direction in the moon's wake for typical interplanetary conditions. 
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S=7 and 0=135°. It should be noted that throughout this paper, angles 
are measured from the X-axis, the direction of the relative velocity of 
the solar wind with respect to the moon. In general the X-direction is 
slightly deflected from the sun-moon line; the deflection angle is 
determined by the following two factors: (1) The direction of the solar 
wind is not strictly along the heliocentric radial direction; the 
deviation from radial direction may be as large as 10°, the average 
direction of the solar wind 15 ' 16 is from east of the sun by 1.5°. (2) 

The moon is orbiting around the sun with a velocity of 30 + 1 km/sec, 
and must be included in calculating the relative velocity of the solar 
wind with respect to the moon. This latter effect is the aberration 
of the solar wind flow which is generally between 3-6°. 

The direction of the wake is "steered* 1 by both the solar wind 
velocity and the imbedded field lines. When 0 = 90° or 180°, the flow 
is symmetrical about the XZ-plane , and the direction of the wake coin- 
cides with that of the upstream solar wind velocity with respect to the 
moon. From Figs. 4 and 5 we can see that when 0 is between 90° and 
180° the wake is deflected slightly toward the direction of the imbedded 
field lines. As the speed ratio S increases the steering of the wake 
direction becomes more dominantly controlled by the solar wind velocity. 

A perspective view of the moon 1 s wake is plotted in Fig. 7 to 
show the ion flux distribution over cross-sections X = 3,6,9,12 and 15. 

We can see that the flow is symmetrical about the XY-plane which is 
parallel to both the solar wind ve locity , JJ Q , and the interplanetary 
magnetic field . The thickness of the wake, measured perpendicular 
to the plane of symmetry, remains constant at one lunar diameter every- 
where. The width of the wake measured in the plane of symmetry increases 
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as its distance from the moon increases. This shows the important role 
the interplanetary magnetic field plays in affecting the flow conditions 
in this wake region. In the absence of the magnetic field, we would see 
a wake flow symmetrical about the X-axis downstream of a sphere rather 
than a deflected wake symmetrical about the XY-plane. 

The ion flux is of primary importance, because this is the quan- 
tity that can be directly measured from space satellites. 9 The distri- 
bution of ion density is useful in computing the electron flow and the 
magnetic field in the moon’s vicinity. Since the direction of the 
forces between charged particles is such as to try to restore neutrality, 
it is characteristic of plasma to try to remain as electrically neutral 
as possible. The ion density distribution acts as a background for the 
electron flow in the moon's wake, and the electron density can deviate 
from the ion density only very slightly. 

Finally we shall discuss the distribution of ion density in the 
vicinity of the moon. Figs. 8 and 9 show the ion density distribution 
compared with the ion flux distribution. The constant density contours 
are quite close to the constant flux contours, and thus the general 
feature of the moon's wake does not change whether we are dealing with 
the density or the flux. 

As to the length of the wake, this author does not attempt to de- 
fine it precisely. A typical variation of ion density and ion flux 
along the X-axis is given in Table 1. From this table it may be stated 
that the length of the moon's wake is of the order of 100 lunar radii. 



12 


Summary 

As the solar wind interacts with the moon, the motion of ion 
particles in the vicinity of the moon is considered as a free molecule 
flow of guiding-center particles. Analytical results are obtained to 
describe the distribution of the ion density and the ion flux in the 
vicinity of the moon. 

The interplanetary magnetic field has a strong effect on the flow 
of the magnetized solar wind plasma. The flow condition is symmetrical 
about a plane passing through the center of the moon and parallel to 
both the solar wind velocity and the imbedded field lines. A sizable 
empty cavity is produced in the near-downstream vicinity of the moon. 

The disturbed flow forms a long slightly curved and deflected wake on 
the darkside of the moon. The thickness of the wake, measured perpen- 
dicular to the plane of symmetry, remains constant at one lunar diameter. 
The width of the wake measured in the plane of symmetry, increases with 
its distance from the moon. The direction of the wake essentially 
follows that of the solar wind velocity with respect to the moon; it 
is also steered by the field lines to bend slightly toward the direction 
of the field lines. As the speed ratio increases the deflection angle 
decreases, and the steering of the wake direction becomes more dominantly 
controlled by the solar wind velocity. The length of the wake is of 


the order of 100 lunar radii. 
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TABLE 1 


Variation of ion density and ion flux along X-Axis (for 
S = 5 and 0 = 135°) 


X 

n/no 

nu/rioUo 

25.0 

0.68849 

0.68801 

50.0 

0.84139 

0.84133 

75.0 

0.89390 

0.89388 

100.0 

0.92033 

0.92032 

125.0 

0.93623 

0.93623 

150.0 

0.94684 

0.94684 

175.0 

0.95443 

0.95443 

200.0 

0.96012 

0.96012 

225.0 

0.96455 

0.96455 

250.0 

0.96809 

0.96809 

275.0 

0.97099 

0.97099 

300.0 

0.97341 

0.97341 

325.0 

0.97545 

0.97545 

350.0 

0.97721 

0.97721 

375.0 

0.97873 

0.97873 

400.0 

0.98005 

0.98005 
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Figure 1 
Figure 2 

Figure 3 
Figure 4 

Figure 5 

Figure 6 

Figure 7 
Figure 8 

Figure 9 


The coordinate system. 

The guiding-center velocity vector and the guiding-center 
distribution function. 

Three different regions in the disturbed flow. 

The distribution of ion flux in the XY-coordinate plane 
for varying direction angle 0. 

The distribution of ion flux in the XY-coordinate plane 
for varying speed ratio S. 

The direction of ion flux in the moon’s wake. (The shaded 
areas indicate the magnitude of ion flux). 

A perspective view of the moon's wake (for S=5, 0 = 135°). 
A comparison of the ion density distribution and the ion 
flux distribution. 

A comparison of the ion density distribution and the ion 


flux distribution. 
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